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ABSTRACT 

We explore the role of X-ray photoevaporation in the evolution and dispersal of vis- 
cously evolving T-Tauri discs. We show that the X-ray photoevaporation wind rates 
scale linearly with X-ray luminosity, such that the observed range of X-ray luminosi- 
ties for solar-type T-Tauri stars (10^^-10'^^ erg s""'^) gives rise to vigorous disc winds 
with rates of order 10^^°-10~^ yr~^. These mass-loss rates are comparable to typ- 
ically observed T-Tauri accretion rates, immediately demonstrating the relevance of 
X-ray photoevaporation to disc evolution. We use the wind solutions from radiation- 
hydrodynamic models, coupled to a viscous evolution model to construct a population 
synthesis model so that we may study the physical properties of evolving discs and 
so-called 'transition discs'. Current observations of disc lifetimes and accretion rates 
can be matched by our model assuming a viscosity parameter a = 2.5 x 10^'^. 

Our models confirm that X-rays play a dominant role in the evolution and disper- 
sal of protoplanetary discs giving rise to the observed diverse population of inner hole 
'transition' sources which include those with massive outer discs, those with gas in 
their inner holes and those with detectable accretion signatures. To help understand 
the nature of observed transition discs we present a diagnostic diagram based on ac- 
cretion rates versus inner hole sizes that demonstrate that, contrary to recent claims, 
many of the observed accreting and non accreting transition discs can easily be ex- 
plained by X-ray photoevaporation. However, we draw attention to a smaller but still 
significant, population of strongly accreting (~10~^ M0 yr~^) transition discs with 
large inner holes (> 20 AU) that lie outside the predicted X-ray photoevaporation 
region, suggesting a different origin for their inner holes. 

Finally, we confirm the conjecture of Drake et al. (2009), that accretion is sup- 
pressed by the X-rays through 'photoevaporation starved accretion' and predict this 
effect can give rise to a negative correlation between X-ray luminosity and accretion 
rate, as reported in the Orion data. We also demonstrate that our model can repli- 
cate the observed difference in X-ray properties between accreting and non-accreting 
T-Tauri stars. 
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1 INTRODUCTION 

Protoplanetary discs are a natural outcome of low mass star 
formation, providing a reservoir of material from which the 
star itself continues to accrete and from which planets may 
later form. The accretion history of a newly formed star, the 
evolution of its disc and the formation of a planetary sys- 
tem are all intimately linked and affected by irradiation from 
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the central star. A lot of attention has recently been paid 
to the final dispersal of protoplanetary discs as this sets the 
time-scale over which planets may form. Observationally, the 
study of the dust content of these discs through the analysis 
of their continuum spectral energy distributions (SEDs), in 
the infra-red (IR), has enormously advanced over the past 
few decades. These observations have indicated discs evolve 
in a way that suggest 'standard' viscous evolution (Hart- 
mann et al. 1998) for the first few Myrs (e.g. Haisch et al. 
2001b), but then rapidly evolve from a disc-bearing (primor- 



© 2002 RAS 



2 J. E. Owen, B. Ercolano & C. J. Clarke 



dial) to disc-less status (e.g. Luhman et al. 2010; MuzeroUe 
et al. 2010). This rapid clearing of the inner disc indicated 
from IR observations has been supplemented with comple- 
mentary observations in the sub-millimetre that show most 
non-accreting stars (WTTs) are devoid of emission out to 
500 AU (Duvert et al. 2000). This suggests the removal 
of material material close to the star (< 1 AU) is correlated 
with the removal of material at larger radius, (> 10 AU); 
and is hence, associated with the entire destruction of the 
protoplanetary disc. 

Furthermore, there have been several observations of 
discs thought to be caught in the act of clearing (e.g Calvet 
et al. 2002; Najita et al. 2007); these 'transition' discs show a 
deficit in emission at near-IR (NIR) wavelengths compared 
to a primordial optically thick disc (one that is optically 
thick all the way to the dust destruction radius), while re- 
turning to the emission levels expected from an optically 
thick disc at longer wavelengths. Currently, there is no clear 
consensus as to what constitutes a transition disc and both 
accreting (e.g. Hughes et al. 2009; Espaillat et al. 2010) and 
non-accreting (e.g. Cieza et al. 2010; Merm et al. 2010) ob- 
jects have been discovered. 

The frequency of these inner-hole sources is approxi- 
mately 10 to 20% of the total disc population for young 
solar-type stars in nearby star forming regions (e.g Strom et 
al. 1989; Skrutskie et al. 1990; Luhman et al. 2010). If one 
assumes that these inner hole objects are indeed 'transition' 
discs (i.e. discs that are caught in the act of dispersal) and 
that the gas and dust evolution go hand-in-hand, then a 10% 
frequency of transition discs suggest a dispersal mechanism 
that operates over a tenth of the optically thick disc lifetime 
(Kenyon & Hartmann 1995). Moreover the colours of such 
discs imply that they clear from the inside out (Ercolano, 
Clarke & HaU, 2010). 

The observations described above have encouraged the- 
orists to propose dispersal mechanism that operate over two 
separate time-scales: a first time-scale of a few million years 
which allows discs to remain optically thick at IR wave- 
lengths, followed by a brief dispersal phase that leaves the 
discs optically thin with colours consistent with those of bare 
photospheres or debris discs. The proposed mechanisms in- 
clude planet formation itself (e.g. Armitage & Hansen 1999), 
grain growth (DuUemond & Dominik 2005), photophoresis 
(Krauss et al 2007), MRI-driven winds (Suzuki & Inutsuka 
2009), and photoevaporation due to extreme ultraviolet ra- 
diation -EUV- (Clarke et al 2001; Alexander et al 2006a,b; 
Richling & Yorke 2000), X-ray radiation (Alexander et al. 
2004; Ercolano et al 2008b, 2009; Gorti & HoUenbach 2009, 
Owen et al. 2010) and far ultraviolet radiation -FUV- (Gorti 
& HoUenbach 2009; Gorti et al. 2009). However there is stiU 
no consensus over which mechanism (or which combination) 
may dominate. 

The recent developments of X-ray photoevaporation 
(XPE) models for T-Tauri discs have yielded encouraging 
results, suggesting that this may be the dominant player 
in disc dispersal. In particular, we have shown that X-ray 
photoevaporation rates exceed the EUV photoevaporation 
rates by two orders of magnitude (Ercolano et al 2009) and 
can easily reproduce the two time-scale behaviour suggested 
by the observations (Owen et al 2010). This model does not 
suffer from the uncertainties in the heating rates that plague 
the FUV model, due to the unknown abundance of PAHs in 



discs and the very uncertain photoelectric yields (see discus- 
sion in Ercolano & Owen 2010). Furthermore, Ercolano & 
Clarke (2010) discussed the role of metallicity in disc disper- 
sal, and found that photoevaporation predicts shorter disc 
lifetimes at lower metallicity, in contrast to planet forma- 
tion which would predict a very strong negative correlation. 
Yasui et al (2009) have presented observation that favour 
a positive correlation between disc lifetimes and metallicity 
and hence which favour XPE; however, it is too early to draw 
any definitive conclusions about the role metallicity plays 
in disc dispersal. Finally, the slightly blue-shifted forbidden 
line spectrum of discs seen almost face on (e.g. Hartigan et 
al 1995; Pascucci et al 2009) is well reproduced by the XPE 
model for all of the low ionisation and atomic diagnostics 
considered (Ercolano & Owen, 2010), an improvement over 
the EUV photoevaporation model which cannot reproduce 
the OI 6300A line intensities (e.g. Font et al 2004). 

More recent observations have focused on understand- 
ing the nature of some known transition discs (TDs) by in- 
vestigating their gas content and accretion properties. Ev- 
idence for accretion and gas inside the dust inner hole has 
sometimes been used to argue against photoevaporation be- 
ing the agent of inner-hole clearing (e.g. Espaillat et al 2010). 
This is however a misconception, perhaps based on ear- 
lier EUV-based photoevaporation models (Alexander et al. 
2006b). In this work we show that the XPE model of Owen 
et al (2010) clearly predicts a transition phase where there is 
a detectable quantity of accreting gas inside the dust inner 
hole and is indeed consistent with the accretion versus inner 
hole size properties of most observed TDs. 

In this paper we use radiation-hydrodynamics simula- 
tions to develop a X-ray photoevaporation model which is 
then used to construct a synthetic disc population to com- 
pare with available observed disc statistics. We discuss the 
X-ray luminosities of TTs in Section 2, while Section 3 pro- 
vides a description of the X-ray photoevaporation model. In 
Section 4 we discuss a viscously evolving photoevaporating 
disc model and the construction of our disc population. In 
Section 5 we present the results from our population synthe- 
sis model while our final conclusions are given in Section 6. 



2 X-RAY LUMINOSITY FUNCTION 

The X-ray luminosity (Lx), is a crucial input for the X-ray 
photoevaporation model, where we take the X-ray luminos- 
ity to span the range 0.1 — lOkeV. As discussed in Section 3, 
the choice of Lx fully describes the mass-loss rate due to 
photoevaporation in a given system. 

The input X-ray spectrum is identical to the spectrum 
used in previous work (Ercolano et al. 2008b; Ercolano et al. 
2009; Owen et al. 2010; Ercolano & Owen 2010). It is a syn- 
thetic spectrum generated by the plasma code of Kashyap 
& Drake (2000) and fits to Chandra spectra of T Tauri stars 
(e.g. Maggio et al. 2007). It includes an extreme ultraviolet 
component - EUV - (13.6-100 eV), soft X-ray component 
(0.1-1 keV) and a hard X-ray component (> 1 keV). Er- 
colano et al. (2009) studied the effect of attenuating this 
X-ray spectrum with neutral columns. They found that at 
columns of order lO'^^ cm"'^ the soft X-ray component was 
screened out and the photoevaporation rates dropped signif- 
icantly. However, for columns ^ 10^° cm~^ the photoevapo- 
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ration rates remained unaffected. Furthermore, Ercolano & 
Owen (20f 0) showed that the X-ray wind itself is optically 
thick to EUV photons: thus when considering photoevapo- 
ration it is only the strength of the soft X-rays that is rele- 
vant. We assume that the soft X-rays are able to reach the 
disc surface as is consistent with the explanation of the OI 
6300A line emission presented in Ercolano & Owen (2010). 
Therefore, we adopt the unattenuated spectrum shown in 
Owen et al. (2010) and expect our results to generalise to 
cases of moderate screening. We will discuss the implications 
of this assumption further in Section 5. 

It has been known for some time that there is a large 
scatter in the X-ray luminosity of T-Tauri stars for a given 
stellar mass (and bolometric luminosity). Therefore, we have 
used the data for the Taurus cluster (Giidel et al. 2007) and 
the Orion Nebular cluster (Preibisch et al. 2005) to build cu- 
mulative X-ray luminosity functions for the Orion sample 
functions (XLFs) for all pre-main sequence stars (including 
both CTTs and WTTs) in the 0.5 to 1 Mq range. The data 
for these two regions shows good agreement at low luminos- 
ity, but they differ at high luminosities, in the sense that the 
Orion data contains a higher fraction of sources with high 
Lx- The difference is most probably not intrinsic to the X- 
ray properties of these two regions, but due to the different 
treatment of strong flares in the two samples. Strong flares 
were excluded in the Taurus data, in contrast to the Orion 
sample where luminosities were averaged over the whole ob- 
servational period. Strong flares are relatively rare and Al- 
bacete Colombo et al. (2007) found that non-flaring sources 
have a median KT — 2.1 ± O.SKeV, compared to flaring 
sources with KT — 3.8KeV. Therefore, due to their higher 
X-ray temperature flares tend to emit most of their radiation 
in the hard X-ray region (as shown by observations of ob- 
jects in a 'flare' state compared to in a 'quiescent' state, e.g. 
Imanishi et al. 2001), where the thermal impact is low, as 
discussed above. For this reason, the Taurus sample is the 
most appropriate for use in the photoevaporation model, 
since it should better approximate the quiescent and there- 
fore softer X-ray luminosity, which as discussed above is the 
relevant input for calculations of the photoevaporation rate. 
Therefore, we adopt the Taurus XLF for the remainder of 
this work. 

We assume that the X-ray luminosity function remains 
invariant throughout the stars' pre-main sequence evolution. 
While it is known that the median of the stellar X-ray lu- 
minosity function does in fact decrease with age due to stel- 
lar spin down (e.g. Hemplemann et al. 1995; Giidel et al. 
1997; Giidel 2004), the evolution of the X-ray luminosity 
during the disc dispersal phase is much smaller (see Fig- 
ure 41, Giidel 2004). Certainly, any evolution of the X-ray 
luminosity for ages up to several tens of Myrs is smaller 
than the observed spread in X-ray luminosities at ~ IMyr. 
Furthermore, there is observational evidence that the X-ray 
luminosities of CTTs and WTTs are significantly different, 
with WTTs in general being more luminous than CTTs (e.g. 
Neuhauser et al 1995; Stelzer & Neuhauser 2001; Flaccomio 



Where the X-ray luminosity is taken to be in the range 0.3- 
lOkeV for the Taurus sample and 0.5-8keV for the Orion sample, 
consistent with our definition which includes both the soft and 
hard component. 
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Figure 1. Cumulative X-ray luminosity functions for the Taurus 
(black soUd Une) and Orion (red dashed Une) clusters for solar 
type stars with mass in the range [0.5,1.0] Mq, including both 
WTTs and CTTs. 

et al. 2003; Preibisch et al. 2005). This has lead to a dis- 
cussion of X-ray emission being 'disturbed' by accretion, 
in terms of either X-ray absorption in accretion columns 
(Gregory et al. 2007) or confinement of the X-ray producing 
corona in accreting systems (Preibisch et al. 2005). Recently, 
Drake et al. (2009) suggested that X-rays may modulate 
accretion through photoevaporation (something they called 
'photoevaporation starved accretion'). Such a scenario may 
be able to account for the difference in the observed X-ray 
luminosities of CTTs and WTTs since more luminous stars 
will lose their disc's first. 

In order to assess the effect of 'photoevaporation starved 
accretion' in explaining the X-ray observations we adopt 
here the null stance; the X-ray luminosity of an individual 
young stellar object (YSO) remains constant in time as our 
models evolve from CTTs to WTTs (i.e. we assume that the 
X-ray luminosity function is fixed). 



3 X-RAY PHOTOEVAPORATION 

While the first self-consistent numerical simulations per- 
formed by Owen et al. (2010) were a significant step for- 
ward in understanding X-ray photoevaporation, the models 
were calculated for only one value of the X-ray luminos- 
ity (2 X lO'^'erg s"'^). A variety of underlying disc models 
were considered, which showed that X-ray photoevapora- 
tion is fairly insensitive to the details of the X-ray 'dark' 
region, due to the fact that the sonic surface is located at 
least several fiow scale heights from the X-ray 'dark'/'bright' 
transition. However, the large observed spread in X-ray lu- 
minosities means that the dependence of photoevaporation 
on the X-ray luminosity needs to be considered. 

Before we discuss the results of a detailed numerical in- 
vestigation of parameter space, we can use simple fiuid me- 
chanics to predict the variation of photoevaporation rates 
with X-ray luminosity. Namely, any axi-symmetric inviscid 
steady state fiow must satisfy the following conditions: con- 
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servation of mass; conservation of specific angular momen- 
tum; hydrostatic equilibrium perpendicular to the stream- 
lines; and the conservation of a Bernoulli type potential 
along each streamline. As discussed in Ercolano & Owen 
(2010) the velocity and temperature in a thermally driven 
wind is intimately linked the value of the effective potential, 
and unlikely to be greatly affected by changes in the X-ray 
luminosity. Furthermore, as shown in Owen et al. (2010) the 
gas temperature can be roughly described in terms of the 
ionization parameter (^ = Lx/nr"^) alone. Thus if we con- 
sider any X-ray heated wind (primordial or a disc with an 
inner hole) and vary the X-ray luminosity, spatially fixing 
the temperature and velocity requires that the ionization 
parameter remains constant. This means that the density 
will scale globally in the flow as n oc Lx- It is then a trivial 
matter to show that the new wind solution we have con- 
structed, with a different X-ray luminosity will still satisfy 
the required conditions described previously. Hence, we can 
use the conservation of mass to show that the X-ray photo- 
evaporation model predicts: 

Mw oc Lx (1) 

for all X-ray flows including both primordial discs and disc 
with inner holes. Also, since the only change is to the global 
density, the normalised mass-loss profiles should also remain 
approximately independent of the X-ray luminosity. Fur- 
thermore, Owen et al. (2010) also showed that there was 
very little variation in mass-loss rate with inner hole radius. 

3.1 Radiation-Hydrodynamic Models 

We have performed radiation-hydrodynamic models to 
study the form of the mass-loss, under the assumption that 
it is driven by X-ray photoevaporation (XPE). We consider 
XPE from discs around solar-type stars for both primordial 
(optically thick discs extending into the dust destruction 
radius) and discs with inner-holes (in both gas and dust). 
We run primordial disc models at six equal logarithmic in- 
tervals between logLx = [28.3,30.8]. We further calculate 
four inner hole models at three different X-ray luminosities 
log Lx = 28.8, 29.8, 30.8 (in addition to the eight models 
already computed by Owen et al. 2010 at logLx = 30.3), 
at radii of ~ 5, 10, 20, 70 AU. We start from an initial den- 
sity structure of a protoplanetary disc surrounding a O.TM© 
star with T'eff=4000 K and Ri,=2.^Rq, taken from the set 
of D'Alessio et al. (2004). We adopt the following elemental 
abundances, given as number densities with respect to hy- 
drogen: He/H= 0.1, C/H= 1.4 x IQ-*, N/H= 8.32 x 10"^ 
0/H= 3.2 X 10"*, Ne/H= 1.2 x 10"*, Mg/H= 1.1 x 10-^ 
Si/H= 1.7 x 10"^ S/H= 2.8 x 10~^ These are solar abun- 
dances (Asplund, Grevesse & Sauval, 2005) depleted accord- 
ing to Savage & Sembach (1996). We use the zeus2D code 

^ While the conservation of the Bernoulli potential does not for- 
mally exist for non-barotropic flows, it is easy to verify that for the 
X-ray winds studied here (i.e. unconfined and thermaUy driven), 
one can find a suitable, well defined (single valued) function P(p) 
along each streamline that allows us to construct an energy in- 
tegral equivalent to the Bernoulli potential for a baratropic flow 
with the identical P(p). By uniqueness, this Bernoulli-type en- 
ergy integral must be a conserved quantity along the streamline 
in the X-ray heated flow. 



(Stone & Norman 1992) to calculate the hydrodynamical 
evolution of the disc, where the gas temperature due to X- 
ray irradiation is parametrised as a function of ionisation 
parameter using the 3D radiative transfer code MOCASSIN 
(Ercolano et al 2003, 2005, 2008). The temperature of gas 
that is not heated by the X-rays is fixed to the dust tem- 
perature and the transition point between the X-ray bright 
and the X-ray dark region occurs at a column of lO'^^ cm~^. 
Such an approximation obviously results in a temperature 
and density discontinuity at this point; however, since this 
transition occurs well below the sonic surface (where the 
mass-loss rate is effectively set) it does not affect the re- 
sulting mass-loss rate. The distribution evolves to a steady 
state, with a bound X-ray dark disc and a thermally driven. 
X-ray bright, transonic photoevaporative wind. The numer- 
ical methods employed here and briefly described above are 
similar to those of Owen et al. (2010) and we refer the reader 
to that article for a more detailed description of the model 
setup. 

3. LI Wmd Rates & Streamline Topology 

The numerical models described above were used to cali- 
brate the analytical relations derived at the beginning of 
Section 3, which we can use to build a synthetic disc popu- 
lation. The main results are summarised as follows: 

(i) As expected, for a star of a given mass the total mass 
loss rate scales almost linearly with the X-ray luminosity. 
This result applies to both primordial and inner hole sources. 
Specifically, 

M. = 6.4xlO-^A(^^^3i|-,y"Moyr-i (2) 

where j4 is a constant taking the value unity for primordial 
discs and 0.75 for discs with inner holes. Figure [2] shows 
the mass loss rates calculated by our grid of radiation- 
hydrodynamics models for primordial discs. Equation [2] is 
also a good fit for inner-hole sources, irrespective of hole 
size, for the cases studied here (inner hole radii between 5 
and 70 AU). Figure|3]shows the results for inner hole sources, 
where the left and right panels show the mass-loss rate as a 
function of hole radius and X-ray luminosity respectively. 

(ii) The cumulative mass-loss rates (i.e. the radial profiles 
of the surface mass-loss rates) for different X-ray luminosi- 
ties are fairly self-similar for primordial discs, and also for 
inner hole sources when one considers the profile as a func- 
tion of {R — Rin), where Ri„ is the inner hole radiu^. Fig- 
ure 13] (left panel) shows the normalised cumulative surface 
mass-loss rates as a function of radius for different X-ray 
luminosities for our primordial disc models. It is clear from 
the figure that higher X-ray luminosities produce a broader 
mass-loss profile; however, this has negligible effect on the 
global disc evolution and we adopt a mean profile, shown 
by the solid black line in the figure. Figure |3] (right panel) 
shows the normalised cumulative surface mass loss rates as 
a function of R-Rin for the inner hole models. The mean 
profile adopted for the inner hole models is shown again as 
the solid black line. 

^ Throughout this work we use R, cf>, z to refer to cylindrical co- 
ordinates and r, 6, ip to refer to spherical polar co-ordinates. 
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Figure 3. Mass-loss rate as a function of inner hole radius (left panel), the points show the results from the simulations (upright triangles 
- Lx = 30.8; circles - results taken from Owen et al. 2010 for Lx = 30.3; squares - Lx = 29.8; downward triangles - Lx = 28.8) and 
the dashed line shows the fit used throughout this work. The right panel shows the mass-loss rate as a function of X-ray luminosity (• - 
~ 5AU; X - Si lOAU; -(- - ~ 20AU; * - 70AU) with the dashed line showing the fit used in the viscous evolution. 
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Figure 4. Left: radial mass-loss profiles for primordial discs. The points represent the results from individual hydrodynamic solutions. 
Right: mass-loss profiles for discs with inner holes plotted as a function of {R— Ri„). The points represent a subset of the 21 hydrodynamic 
models calculated. All models show the same general profile with higher luminosities and smaller inner holes giving rise to slightly broader 
profiles. In both panels the solid black line represent the fits used to both primordial and inner hole discs. 



From the results presented in tills section It Is Immedi- 
ately apparent that the mass-loss properties for a star of a 
given mass are completely controlled by the X-ray luminos- 
ity. We can therefore construct a population synthesis model 
that Includes only viscosity and XPE (with appropriate Ini- 
tial conditions, see Section 4.1). 



4 PHOTOEVAPORATING VISCOUS DISCS 

The evolution of the surface density of a photoevaporatlng 
and viscously evolving disc can be described In one dimen- 
sion using the formalism of Lynden-Bell & Prlngle (1974): 



as 

dt 



3__d_ 
RdR 



tn,{R,t) (3) 
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X-ray Luminosity [erg s ^] 

Figure 2. Total mass-loss rates as a function of X-ray luminosity. 
Black squares represent the individual hydrodynamic models, the 
solid line represents the fit used in this work (Equation (JJ, the 
dashed line represents the linear analytical prediction. 



where !^(R) describes the viscosity term and T,w{R,t) rep- 
resents the mass-loss due to photoevaporation calculated in 
the previous section. Before moving on to discuss our choice 
of initial conditions and viscosity law, it is worthwhile dis- 
cussing the qualitative evolution of a photoevaporating vis- 
cously evolving disc. Equation [3] is perhaps easiest to un- 
derstand in its integral form (i.e dR2iiR x Equation [3]) , 
which tells us: 



dt 



(4) 



where Md is the disc mass and A/* is the accretion rate onto 
the star due to viscous transport. Equation [4] tells us im- 
mediately that there are two phases of disc evolution: when 
A/« > Mji, the disc evolution is dominated by the viscous 
transport of angular momentum and hence associated ac- 
cretion rather than the removal of gas through photoevap- 
oration. During this stage the disc will behave in a similar 
manner to a standard viscously evolving disc, without pho- 
toevaportion especially when M, S> Mw- However, when 
M-w > M, it is now photoevaporation that dominates the 
evolution of the disc, opening a gap in the disc and then fi- 
nally removing the remaining disc material until the disc is 
dispersed. The accretion rate of a disc evolves on the disc's 
viscous time ti, at the outer radius, which is of order the 
disc's lifetime (Lynden-Bell & Pringle, 1974; Pringle, 1981; 
Ruden, 1993). The transition from a viscously evolving disc 
to a clearing disc occurs when a gap opens in the inner disc 
and the material inside the gap drains onto the central star. 
Since, the time-scale for the inner disc to drain is the viscous 
time at the point the gap opens (which is typically much less 
than the global viscous time of the disc), the transition from 
a viscously accreting disc (i.e. one that would observation- 
ally be classified as an optically thick accreting primordial 
disc) to one that is being dispersed, occurs on a time-scale 
much shorter than the disc's current lifetime (Clarke et al. 
2001; Ruden 2004). 

Hence, the evolution of a viscously evolving photoevap- 



orating disc can naturally explain the observations of discs 
which appear to evolve under the effect of viscous trans- 
port alone, but are then rapidly dispersed. A key test for 
a model of viscously evolving photoevaporating discs is to 
check whether it can self consistently explain the observa- 
tional statistics for both the evolution of primordial discs and 
those thought to be in the act of clearing (i.e. 'transition' 
discs) and the time-scales associated with their evolution. 
As such, any model of disc evolution depends on the choice 
of viscosity law and the initial conditions associated with 
disc evolution. 



4.1 Viscosity Law and Initial Conditions 

We now turn our attention to the choice of initial conditions 
and viscosity laws that are required to solve Equation[21 We 
adopt the form ^{R) — voR, where uq — acs'^/ft, which we 
always evaluate at lAU. Such a choice is appealing both 
observationally and theoretically, since it predicts a surface 
density scaling as E oc R~^ which is supported by observa- 
tions (e.g. Hartmann et al 1998; Andrews et al. 2010) and is 
consistent with a constant a disc which is mildly flaring i.e. 
H/R oc R^^*. Furthermore, we adopt the zero-time similar- 
ity solution of Lynden-Bell & Pringle (1974) for which the 
initial surface density distribution takes the form: 



E(i?,0) 



2tvRRi 



exp{-R/Ri) 



(5) 



where Md{0) and Ri are the initial disc mass and a scale 
radius describing the exponential taper of the disc's outer 
region . These initial conditions, together with the viscosity 
parameter a and the X-ray luminosity fully determine the 
evolution of the disc through Equation [S] 

At this point it is useful to construct a 'null model' 
of viscous evolution without photoevaporation that when 
combined with the observed X-ray luminosity function can 
explain the observed decline in disc fractions with age. In 
effect we are asking whether there is a universal set of disc 
viscous parameters which can explain the variation in disc 
lifetime from cluster to cluster purely in terms of the ob- 
served spread in X-ray luminosity. As discussed above the 
basic principle of all photoevaporation models is that discs 
should evolve viscously, hardly noticing the effects of photo- 
evaporation, until the mass accretion rates in the discs have 
fallen to a value that is comparable to the photoevaporation 
rat(EI, at which point the remaining disc material is rapidly 
cleared. 

Viscous evolution alone with v cc R predicts that accre- 
tion rates evolve as: 



M,{t) = M.(0) 1 + 



-3/2 



(6) 



where t^, is the viscous timescale at Ri. This evolution 
should therefore be observed in discs before XPE sets in. We 
can equate the fraction of disc-bearing pre-main sequence 
stars at a given time (/d), with the fraction of stars in the 



* Owen et al (2010) showed that XPE can create a gap in a disc 
only when the accretion rates onto the star are approximately an 
order of magnitude lower than the photoevaporation rates. 
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Figure 5. Accretion rates as a function of time for the 'null' disc 
model. The errors in accretion rate are estimated from the \/N 
errors in the disc fractions. The solid line shows a fit of Equation|6] 
used to determine suitable values of Mj{0) and ti^. 



X-ray luminosity function that have luminosities less than a 
cut-off X-ray luminosity Lc{fd)- In order for objects with X- 
ray luminosities equal to Lc to be about to lose their discs, 
the viscous accretion rate must be equal to M^i^Lc) at this 
point. We have performed this simple exercise using the disc 
fractions for nearby clusters compiled by Mamajek (2009), 
the Taurus X-ray luminosity function and the XPE theory 
developed above. The result is shown in Figure[51 where each 
point represents the current accretion rate cut-off in a clus- 
ter implied by XPE. We scale our results on disc fractions 
assuming an initial close binary fraction of 14% by consid- 
ering that the 0.3Myr old cluster NGC 2024 (Haisch et al. 
2001a), which shows a disc fraction of 86%, is too young 
for any disc to have been destroyed by photoevaporation or 
planet formation, but only through binary interactions. The 
solid line in the plot represents a suitable fit of Equation [6] 
to the data, from this fit we can extract an initial accre- 
tion rate of M*(0) = 5 x 1O~*M0 yr~^and a viscous time of 
t,y = 7 X lO^yr. From these two values we can calculate an 
initial disc mass of Md{0) ~ O.OTM©, which is similar to the 
canonical value (10% of the stellar mass) at which viscous 
angular momentum transport takes over from self-gravity. 

Along with giving us appropriate initial conditions for 
our disc population model, the above also provides a strin- 
gent test of the hypothesis that X-rays are key to disc evo- 
lution and dispersal. If the X-rays were not the dominant 
dispersal mechanism, there is no a priori reason to expect a 
'null' model constructed only using knowledge of the X-ray 
luminosity function and observed disc fractions to reproduce 
a plausible evolution of the accretion rates seen in CTTs, 
both in terms of the time exponent (in Equation [6| and 
its initial value. In fact, increasing or decreasing the spread 
about the median of the Taurus X-ray luminosity function 
by a factor of ~5 or greater makes a fit of Equation[6]to the 
'null' model impossible. Although this agreement could be 
fortuitous, it is reassuring that the X-ray luminosity func- 
tion, disc lifetimes and accretion histories are consistent with 
our XPE hypothesis. 

Furthermore, in order to uniquely specify the viscous 
evolution we must pick suitable values of a and Ri which in 
turn specifies pq. While any combination of Ri and a that 
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Figure 6. Primordial disc fraction as a function of time (solid 
black line) from our XPE population calculated using 500 disc 
models with the mass-loss rate determined by randomly sam- 
pling the Taurus X-ray luminosity function, all discs evolve from 
a single set of initial conditions. The points are observed disc 
fractions compiled by Mamajek (2009). Model disc fractions have 
been scaled to account for disc destruction by close binary in- 
teraction. Note that the structure of the black line reflects the 
structure in the X-ray luminosity function. 



give the required viscous time will reproduce the same 'null' 
viscous model (i.e. non-photoevaporating) , a disc evolution 
model that includes photoevaporation mildly breaks this de- 
generacy. By performing a fit (by eye) of a viscously evolving 
photoevaporating disc population to the disc fractions used 
to derive the 'null' model we obtain values of a = 2.5 x 10"'^ 
and Ri = 18AU although due to the large scatter in the disc 
fractions and only mild breaking of the a, Ri degeneracy, 
we estimate errors on these values of a factor of 2-3. In Fig- 
ure[6]we show the obtained fit to the disc fractions compiled 
by Mamajek (2009) for our viscously evolving photoevap- 
orating disc model with a single set of initial conditions: 
Md{0) = O.O7M0, a = 2.5 X 10"^ and Ri = 18AU. 

Finally, Figure [7] shows the predicted evolution of ac- 
cretion rate with time for individual disc models (each line 
shows the evolution of accretion rate with time for one disc 
model, for a given value of the X-ray photoevaporation rate). 
This clearly shows that a spread in accretion rates at late 
times > IMyr does not necessary require a spread in initial 
conditions as some authors have suggested is necessary out- 
side the XPE framework (e.g. Armitage et al. 2003; Alexan- 
der & Armitage, 2009). Observationally there is however, a 
spread in accretion rates seen at early times <lMyr (e.g. 
Hartmann et al. 1998). This variability cannot be fit with 
a single set of initial conditions since photoevaporation has 
had no time to act on the disc. The failure to match the 
spread in accretion rates at early times < IMyr may not 
be surprising, since at early times the angular momentum 
transport mechanism may be dominated by self-gravity and 
accretion may be episodic (e.g. Lodato & Rice, 2005). In 
reality the zero time point in our models corresponds to 
the point at which the transition from a self-gravitating to 
a viscous disc occurs, as indicated by the determination of 
Afd(O) = O.IM*. This may explain, why our initial accretion 
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Figure 7. Each solid line represents the evolution of accretion 
rate as a function of time for an individual disc model under- 
going X-ray photoevaporation. All models have the same initial 
conditions, the X-ray luminosity is the only variable. 



rate is significantly lower than some of the accretion rates 
msEisured for the youngest objects (Hartmann et al. 1998). 

Furthermore, better agreement for the early time data 
could easily be obtained by assuming a range of initial sur- 
face density profiles. Indeed our choice of a self-similar sur- 
face density distribution at zero time has no physical ba- 
sis other than convenience (in fact it would be extremely 
surprising for discs to be born with the surface density dis- 
tribution of the zero time similarity solution of Equation [3] 
given viscous angular momentum transport is not a key pro- 
cess during the disc formation stage). We emphasise, that 
any initial surface density distribution with the same initial 
viscous parameters Md(0), tv and a will tend to the same 
evolution after a few viscous times as shown in Lynden-Bell 
& Pringle (1974). For this reason, we do not attempt to ex- 
plain disc evolution at early times here, since we are mainly 
interested in the question of disc dispersal, for which only 
the viscous evolution phase at > IMyr is relevant. 



5 RESULTS AND DISCUSSION 

We have used the methods and initial conditions derived 
in the previous section to construct a population synthesis 
model for the evolution of discs dominated by viscosity and 
XPE. We have calculated 500 disc models based on a ran- 
dom sampling of the Taurus X-ray luminosity function. Our 
disc evolution models are computed by solving Equation [3] 
numerically using the method set out in Owen et al. (2010), 
following the evolution of the disc until the disc is cleared 
to lOOAU. At radii larger than lOOAU the photoevaporation 
rates are extremely uncertain and we cannot be confident in 
results that continue the evolution beyond this radius. How- 
ever, for the sake of completeness we will discuss the possible 
qualitative evolution of these remnant discs in Section 5.3. 
We now turn our attention to some specific predictions from 



these models and, where possible, compare them with ob- 
servations. 

5.1 Photoevaporation starved accretion 

As discussed in Section 2, Drake et al. (2009) suggested that 
coronal X-rays suppress the accretion flow onto young solar- 
type stars through the driving of a photoevaporating wind. 
This photoevaporation starved accretion phase can explain 
the tentative negative correlation between mass accretion 
rate and stellar X-ray luminosity reported by Drake et al 
(2009). Moreover the reduction in disc lifetime in strong X- 
ray sources can explain the observation that the X-ray lumi- 
nosities of accreting T Tauri stars are systematically lower 
than those of non-accretors. 

The XPE models of Owen et al (2010) clearly show 
that there is indeed a phase in the disc evolution (before the 
opening of the gap) where the effects of this 'starving' are 
apparent in the radial dependence of the accretion rate. In 
Figure [8] we compare the accretion rate and surface density 
profiles of the median disc model undergoing XPE, 0.5Myr 
before the gap opens, against those of a disc which is only 
subject to viscous evolution. Inside 70 AU the accretion rate 
drops before it reaches the star compared to the standard 
case where the accretion rate tends to a constant throughout 
the entire disc. This can be compared to the EUV photoe- 
vaporation model: in this case the mass-loss profile is nar- 
rowly peaked between 1 — lOAU (for solar type stars) and the 
total mass-loss rate is considerably less (~ lO^^'^M© yr^""^). 
This results in a shorter and much less pronounced period 
of 'starving' i.e. the disc is only affected inside a few AU. 
In contrast, the photoevaporation starved accretion lasts for 
~20-30% of the disc lifetime in the X-ray model, with signif- 
icant consequences for global disc evolution: i.e. a flattening 
of the surface density profile and a significant drop in the 
accretion rate through the disc. 

In Figure [5] we show the evolution of the surface density 
for the median disc model (i.e. the disc with the median X- 
ray luminosity of 1.1 x 10"^'' erg s~^ which corresponds to a 
photoevaporation rate of 7.1 x 10~^ Mq yr~^) undergoing 
the stages of gap-opening and final clearing. This shows the 
drop in surface density through the disc between 1 — 70AU 
before the gap opens (as shown in Figure Moreover the 
broad photoevaporation profile also causes the steady ero- 
sion of the disc during the draining of the inner hole, so that 
the hole (though opening at ~ 3AU) roughly doubles in size 
during inner hole draining. Once the inner disc has com- 
pletely drained there is a rapid clearing of the disc out to 
10 — 20AU, because this region was previously depleted dur- 
ing the photoevaporation starved accretion phase. The fast 
clearing phase slows down once the inner hole reaches radii 
less affected by photoevaporation starved accretion. This re- 
sults in a a total clearing time of roughly 10 — 20% of the 
disc lifetime, which is consistent with transition disc statis- 
tics. This can be compared to the EUV models of Alexander 
et al. (2006b), which resulted in a transition phase approx- 
imately 3% of the disc's lifetime, and those of Clarke et al. 
(2001) which clear the outer disc on time-scales of the order 
of the disc lifetime. 

The negative M-Lx correlation reported by Drake et 
al (2009) is a simple consequence of the fact that discs 
with higher X-ray luminosities produce more vigorous winds 
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Figure 8. Surface density (top) and accretion rate (bottom) ra- 
dial profiles. The solid line represents a disc undergoing no photo- 
evaporation and the dashed line represents the median disc model 
0.5Myr before a gap opens in the disc due to XPE. 
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Figure 11. Time evolution of the median X-ray luminosity in 
the synthetic cluster for CTTs (solid line) and WTTs (dashed 
line). The dotted line represents the transition X-ray luminosity 
between a CTTs and a WTTs as a function of time. The data 
points are taken from the compilation of Giidel et al. (2007) in 
the Taurus cluster, filled circles refer to CTTs and filled squares 
refer to WTTs. 
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Figure 9. The evolution of the disc's surface density during the 
disc clearing phase. The first line shows the zero time surface 
density profile, the next shows the profile at 75% of the discs 
lifetime (~ 3.5Myr) and the remaining lines show the surface 
density at 1% steps in disc lifetime. 



causing the disc's accretion rates to be lower than those for 
discs with a less vigorous wind. Thus a negative M-Lx cor- 
relation is expected for clusters with a relatively narrow age 
range. This effect is however counteracted by the fact that 
discs with lower X-ray luminosities take longer to evolve, 
spending more time at lower accretion rates compared to 
high luminosity objects. Therefore, if the X-ray luminosity 
was compared to accretion rates for an entire disc popu- 



lation's lifetime, in clusters with very large age spreads a 
positive correlation might be expected. 

In Figure [10] we show plots of M-Lx for our synthetic 
'cluster' members selected to give clusters with various nar- 
row age ranges (upper panels and lower left panel) along 
with the a plot of the total disc population (bottom right 
panel). As expected, the XPE model predicts a positive M- 
Lx correlation for clusters with large age spread and a neg- 
ative correlation for clusters with a narrow age spread. As 
a comparison, the age spread in Orion is roughly 2-3Myr 
(Haisch et al. 2001b), which explains the observation of a 
negative correlation by Drake et al. (2009). 

The observation of systematically higher X-ray lumi- 
nosities of non-accreting TTs (WTTs) compared to ac- 
creting TTs (CTTs) can also be explained by our models 
in terms of photoevaporation starved accretion. Figure 1111 
shows the time evolution of the median X-ray luminosity 
of the CTTs (solid line) and WTTs (dashed hne) popula- 
tions for our model compared to the data compiled by Giidel 
et al. (2007) for the Taurus cluster (black circles and red 
squares). We note the good agreement between the model 
predictions and observations. The dotted line represents the 
critical X-ray luminosity as a function of time which in our 
model separates CTTs and WTTs and corresponds to the 
X-ray luminosity of objects which have just opened a gap 
and have begun the clearing phase. The number of anoma- 
lous objects (i.e. CTTs above the line and WTTs below the 
line) is small at ages 1 Myr. Given uncertainties in age 
determinations (particularly at < a Myr) the agreement be- 
tween the observations and predictions is very encouraging. 
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Figure 10. Synthetic observations of X-ray luminosity and accretion rates of clusters with a uniform age spread, the synthetic observations 
are shown as points while a linear fit to the data is shown as the solid line. The top left and right panels and the bottom left panels 
represent young clusters with uniform age spreads of 0.5-3.5Myr, 1.5-4.5Myr and 4-8Myr, respectively. The bottom right panel represents 
the disc population observed over the entire evolution. 



5.2 The nature of transition discs: accreting and 
non-accreting 

The recent observations of a class of transition (inner hole) 
discs with residual gas inside the inner dust radius and 
with signature of accretion has prompted some authors to 
question the viability of photoevaporation as the formation 
mechanism for the inner hole in some sources (e.g. Cieza 
et al 2010). Previous EUV-driven photoevaporation mod- 
els (e.g. Alexander & Armitage 2009) indeed predicted that 
at the time of gap opening the surface density of the gas 
in the inner disc and the accretion rates due to the inner 
disc draining onto the star should be undetectable. This 
is however not the case for XPE, mainly due to the fact 
that the virind rates can be two orders of magnitude higher 
than the EUV-driven rates, meaning that at the time of gap 
opening the mass of the draining inner disc and the accre- 
tion rate onto the star of the inner disc material remain 
detectable for a non- negligible amount of time. For the disc 
population generated in this work. Figure 1121 shows that 
the accreting inner holes and non-accreting inner-holes {M 
< 1 X 1O~^^M0 yr~^) are in general equally likely out to a 



radius of Rin ~ 5AU. Clearly as the transition disc is further 
photoevaporated and its inner radius moves out the accre- 
tion signatures onto the star become less evident and non- 
accreting inner holes dominate at radii larger than 20 AU. 
The total integrated ratio out to lOAU (the radius probed 
by 24/im emission around solar-type stars) is found to be 
25% accreting and 75% non-accreting for the entire popula- 
tion. We caution that this is not equivalent to the observed 
fraction of accreting to non-accreting objects in a individ- 
ual cluster, where the cluster age should also be accounted 
for. In young clusters the transition disc population is dom- 
inated by high X-ray luminosity objects which give rise to a 
considerably longer accreting inner hole phase. In contrast 
this ratio is much lower in old clusters where the transi- 
tion disc population is dominated by low X-ray luminosity 
objects that have very short accreting inner-hole phases. 

It is perhaps worth noting at this point that the detec- 
tion of a 'transition' disc observationally is made through 
observation of the dust continuum spectral energy distri- 
bution (SED). Alexander & Armitage (2007) examined the 
behaviour of dust in a photoevaporating disc, finding that, 
under the action of dust drag, the time-scale for dust grains 
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Figure 12. Probability distribution function of accreting (solid 
line) and non-accreting (dashed line) transition discs with inner 
hole radius, they are both scaled so that the integrated sum of 
both distributions is 100%. 



to drain onto the star is of order lO'^yrs, after the gap opens, 
approximately two order of magnitudes faster than the gas 
draining time-scales. This means that an observer would cer- 
tainly see a significant drop in opacity in the inner disc im- 
mediately after a gap has opened, while the gas will still 
linger in the inner dust disc for the duration of its viscous 
draining time-scale of ~ lO^yrs. 

We have used our population synthesis model to inves- 
tigate the accretion rate versus inner hole size evolution for 
transition discs created by XPE under the assumption of im- 
mediate dust clearing at the time of gap opening. Figure [131 
shows the probability distribution of the disc models in the 
M-Rin plane. The symbols represent a sample of observa- 
tions of solar-type objects classified as 'transition' discs by 
Espaillat et al. (2007a,b, 2008,2010) - Red Circles, Hughes et 
al. (2009,2010) - Red Squares, Kim et al. (2009) - Red Dia- 
monds, Calvet et al. (2005) - Black Diamonds, Merm et al. 
(2010) - Black Squares and Cieza et al. (2010) - Black Trian- 
gles (Although Cieza et al. 2010 do not fit for the inner-hole 
radius they list as transitional sources those discs that have 
a deficit of emission in the Spitzer IRAC bands; therefore 
we conservatively estimate an inner hole radius of < lOAU 
for all their sources). 

It is immediately clear from the figure that there is a 
population of large inner hole, strongly accreting transition 
discs that cannot have been created by XPE. Gap open- 
ing by a giant planet or grain growth is perhaps the most 
plausible explanation for these objects. However, there is a 
significant number of discs with inner holes that are consis- 
tent with an XPE origin. Furthermore we note the lack of 
observations of non-accreting 'transition' discs with holes at 
radii greater than 20AU, where our model predicts a signifi- 
cant populations (although several non-accreting discs with 
large holes have been detected in different mass ranges e.g. 
Merm et al. 2010, where the observations probe different ra- 
dial scales). The observations are still rather sparse and it is 
currently not possible to say whether the observed popula- 
tion of transition discs is a true representation or an artifact 
of observational selection effects. 
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Figure 13. Transition disc probability map in the M-Rin plane. 
We have applied an observational cut-off of 10~^^Mq yr^'^in ac- 
cretion rate and objects below the cut-off are all shown at the 
bottom of the diagram. We show a representative sample of solar- 
type objects classified as 'transition' discs; data from taken from 
Espaillat et al. (2007a,b,2008,2010) - Red Circles, Hughes et al. 
(2009,2010) - Red Squares, Kim et al. (2009) - Red Diamonds, 
Calvet et al. (2005) - Black Diamonds, Merfn et al. (2010) - Black 
Squares and Cieza et al. (2010) - Black Triangles, where error-bars 
are listed they are shown. We also show two model tracks corre- 
sponding to a Idex spread in X-ray luminosities about the median 
model in magenta. 



One obvious consequence of an X-ray photoevapora- 
tion mechanism is that the properties of transition discs 
should be correlated with the X-ray luminosity, something 
no other model of photoevaporation or 'transition' disc ori- 
gin would predict. In Figure [T^ we show two such correla- 
tions namely the inner hole radius (left panel) and accretion 
rate (right panel) against X-ray luminosity, considering only 
accreting 'transition' discs (i.e. those with an accretion rate 
> 1 X 1O~^^M0 yr~^). The plots have been generated by ran- 
domly sampling (in time) the accreting transition phase of 
each disc model several times, and should therefore provide 
a reasonable estimate of both the general form of the corre- 
lation plus the associated scatter. Clearly, since discs with 
higher X-ray luminosities open gaps earlier and at higher ac- 
cretion rates a strong positive correlation between Lx and 
M would be expected and is reproduced by the models as 
shown. 

Furthermore, we also recover a positive correlation be- 
tween inner hole radius and X-ray luminosity for accreting 
transition discs as shown in the figure. We can understand 
this by remembering that while the inner disc is still drain- 
ing the inner hole radius is also being eroded outwards as 
shown in Figure (9) Since, for more X-ray luminous objects 
the region being eroded is more depleted due to photoevapo- 
ration starved accretion and the magnitude of the mass-loss 
is higher, the inner-hole radius is able to be eroded to larger 
radius during the accreting phase. This explains the posi- 
tive slope of the right hand extent of the symbols in the left 
panel of Figure [14l 
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Figure 14. Simulated observations of accreting 'transition' discs. The left hand panel shows X-ray luminosity plotted against inner-hole 
radius. The right hand panel shows X-ray luminosity plotted against accretion rate where an observational cut-off of 10~^^Mq yr~^has 
been used. 



5.2.1 Consequences of a different X-ray luminosity 
function 

The X-ray luminosity function is a crucial input into the 
XPE model. In this work we have used the Taurus X-ray lu- 
minosity function as this best represents the quiescent X-ray 
flux the disc sees throughout its lifetime. However, here we 
discuss the consequences of a different X-ray spectrum that 
may be incident on the disc. If the incident spectrum the disc 
sees is systematically harder than the one used (through at- 
tenuation of the X-ray spectrum by large neutral columns 
close to the star) or softer than the spectrum used, the qual- 
itative behaviour of the disc population would remain the 
same, since one can vary the initial condition to fit. There- 
fore, it is important to assess the qualitative changes relating 
to some of the predictions relating to transition discs. As dis- 
cussed in Section 2 it is only the soft X-rays (0.1-lkeV) that 
have any thermal impact, then the result of a changing spec- 
trum can just be considered to be a change in the soft X-ray 
luminosity incident on the disc. Therefore, overall mass-loss 
rates can simply be scaled to this new harder/softer spec- 
trum. 

In Section 4.1 we argued that any X-ray luminosity 
function with a similar spread can have a 'null' model con- 
structed to fit the observations of disc fraction. Thus, a 
harder or more attenuated spectrum would result in lower- 
mass loss rates, which would require a lower initial accretion 
rate, compared to a softer spectrum, which would require a 



higher initial accretion rate to explain the observed disc frac- 
tions. Provided this effect is systematic across all values of 
Lx (i.e. in doesn't change the spread of X-ray luminosities) 
then the consequences for the predicted 'transition' disc pop- 
ulation can be considered. For the harder/more attenuated 
spectrum a lower initial accretion rate is required implying 
a lower initial disc mass, and therefore a smaller population 
of accreting transition discs, with lower accretion rates and 
smaller inner holes. The converse is true for a softer spec- 
trum which requires a higher initial accretion rate and hence 
larger initial disc mass, giving rise to a larger population of 
accreting transition discs, with higher accretion rates and 
larger holes. 



5.3 Final clearing of the disc 

Our mass loss rates are only accurate out to ~100 AU, since 
at this point other clearing mechanisms (e.g. FUV photoe- 
vaporation) may become dominant. For this reason, we have 
stopped our viscous evolution models once the inner hole 
reaches a size of 100 AU. If we were to extrapolate our mod- 
els to clearing beyond lOOAU using Equation 2, we would 
find that in about 10% of cases (sources irradiated by a 
low X-ray flux) the final clearing time-scales would exceed 
lOMyr, resulting in a population of long-lived discs with 
large inner holes. These discs may not survive long if FUV 
photoevaporation is efficient at large radii as suggested by 
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Gorti & HoUenbach (2009). However, given the uncertain- 
ties in the FUV rates, we also consider the case where these 
objects survive for a long enough time for their dust con- 
tinuum emission to be observed. The spectral energy distri- 
bution (SED) of these cold massive discs should be similar 
to that of young debris discs, where a debris disc model is 
normally considered to be a single temperature belt of op- 
tically thin dust (thought to be constantly replenished by 
collisions between planetesimals) at a given radius from the 
star (Wyatt 2008). This suggests that some of the sources 
that are currently classified as debris discs may in fact be 
'XPE relics'. 

We have used the radiative transfer code of Whitney et 
al (2003a,b) to calculate the SED of a typical XPE relic us- 
ing the standard input disc structure and dust properties of 
the code. The properties of these XPE relics were taken from 
the end point of the lowest X-ray luminosity (and hence the 
longest lived) model with an inner radius of 100 AU, outer e- 
folding radius of 310AU and a mass of 7x lO~^M0in gas (and 
a dust to gas mass ratio of 0.01). While such a disc is likely 
to be settled in the dust, we computed three models: flat 
dust distribution, fully mixed dust distribution and a disc 
where Hdust[R) = Q.lHgas{R)- Figure [TSlshows a plot of the 
fractional excess of the disc compared with the stellar pho- 
tosphere at 24/im and 70/im, where we compare our models 
with observations of objects classified as young debris discs 
around solar- type stars (Wyatt 2008). It is clear from this 
figure that XPE relics with a degree of dust settling share 
the same space in the excess-excess plot as sources currently 
classified as young debris discs. The fully-mixed discs show 
a 70/im excess which is probably too large to be classified 
as a debris disc, although it is extremely unlikely that any 
disc could survive for > lOMyr without undergoing dust 
settling. We also note that the predicted 850/^m emission 
from the evolved XPE relics falls below the current detec- 
tion limits at 50 pc (Andrews & Williams 2005), and thus 
these relics would not contradict previous sub-mm observa- 
tions of WTTs which show that most are devoid of emission 
out to 500 AU (Duvert et al 2000). Only ALMA will be able 
to separate these large massive discs from canonical debris 
discs, and hence confirm or dismiss the existence of this pro- 
posed class of objects, placing constraints on the role of FUV 
photoevaporation at large radius. 



6 CONCLUSIONS 

We have used radiation-hydrodynamic calculations of X-ray 
photoevaporated discs coupled to a viscous evolution model 
to construct a population synthesis model, with which we 
have studied the physical properties of primordial and tran- 
sition discs. The initial conditions and viscosity law are con- 
strained using recent observations of disc fractions in nearby 
clusters (Mamajek 2009). We require a viscosity coefficient 
of approximately a — 2.5 x 10~'^ to match the observational 
constraints. 

Our main conclusions can be briefly summarised as fol- 
lows: 

(i) X-rays play a major role in the evolution and dispersal 
of discs around solar-type stars, driving vigorous photoevap- 
orative winds whose rates scale linearly with the X-ray lu- 



E 



24(im Excess 

Figure 15. Plot of 70/im excess above the photosphere against 
24/im excess for disc models with a large inner hole (100 AU) 
and low X-ray luminosities and therefore a predicted remaining 
lifetime > lOMyr The circles show discs close to edge-on and 
increasing 24/xm excess represents the evolution from a flat to 
fully mixed dust distribution. Similarly the diamonds show the 
disc models close to face-on. Also shown as open squares is ob- 
servations of young solar-type stars (C lOOMyr) with IR excesses 
classified as debris discs, taken from (Wyatt 2008). 



minosity, which are in the range of observed accretion rates 
for T-Tauri stars. 

(ii) We have constructed a 'null' accretion disc model us- 
ing only knowledge of the observed disc fractions and X- 
ray luminosity functions under the assumption that discs 
are dispersed through X-ray photoevaporation. This 'null' 
model shows very good agreement with observed accretion 
rates in YSOs as well as their evolution with time, providing 
further independent confirmation of the viability of X-ray 
photoevaporation as a dominant dispersal mechanism. 

(iii) X-rays suppress accretion by preventing accreting 
material from reaching the star, since this material is re- 
moved through photoevaporation. This 'photoevaporation 
starved accretion' (Drake et al 2009) produces a negative 
correlation between X-ray luminosity and accretion rate for 
clusters with a relatively narrow age spread, in agreement 
with the observational correlation reported by Drake et al. 
(2009) in the Orion data. 

(iv) Our models successfully reproduce the observation 
that WTTs are systematically brighter in X-rays than 
CTTs. Whereas this has previously been interpreted as mod- 
ification of X-ray emission or detectability by the presence 
of discs, our models support Drake et al. (2009) in revers- 
ing the causal link - i.e. that disc lifetimes are regulated 
by XPE and hence discless stars are on average those with 
higher Lx- 

(v) A large fraction (>~ 50%) of observed transition discs 
can be easily explained by X-ray photoevaporation. There is 
however a population of strongly accreting transition discs 
with large inner holes > 20AU that lies outside the M — Ri„ 
region predicted by our models, suggesting that alternative 
mechanisms are responsible for their inner hole (e.g. binary 
interaction, grain-growth and/or planet formation). 



© 2002 RAS, MNRAS OOO.mfTSi 



14 J. E. Owen, B. Ercolano & C. J. Clarke 



(vi) A fraction of currently observed objects classified as 
young debris discs (on the basis of their excesses at 
and 70yum) may in fact be the relics of X-ray photoevapored 
discs, which are predicted to be long-lived (>10 Myr) for low 
X-ray luminosity sources. Future mass determination with 
ALMA are necessary to shed light onto the nature of these 
objects. 
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